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Abstract 

For a class of scalar fields including the massless Klein-Gordon field 
the general relativistic hyperboloidal initial value problems are equivalent 
in a certain sense. By using this equivalence and conformal techniques 
it is proven that the hyperboloidal initial value problem for those scalar 
fields has an unique solution which is weakly asymptotically flat. For data 
sufficiently close to data for flat spacetime there exist a smooth future null 
infinity and a regular future timelike infinity. 



1 Introduction 

1.1 General remarks 

A large open problem of classical general relativity is the characterization of the 
structure of a spacetime by initial data. The flat case, Minkowski spacetime, 
is geodesically complete. To the other extreme the singularity theorems by R. 
Penrose and S. Hawking show that the spacetime cannot be geodesically complete 
if the data are large . 



In the last years there has been remarkable progress in describing what happens if 
one goes from data for flat space to large data: The future of small data evolving 
in accordance with the Einstein equation with various matter models as sources, 
vacuum and Einstein-Maxwell- Yang-Mills, looks like the future of data for flat 



space ||10|, |18l. Nevertheless many problems are still unsolved. 



*This work is too a large extent part of my Ph. D. thesis which has been done at the 
Max-Planck-Institut fiir Astrophysik, Postfach 1523, D-85740 Garching 



Those results were significantly improved by D. Christodoulou for spherically 
symmetric models with a massless Klein-Gordon scalar field as source. He was 
able to relate properties of the initial data to properties of singularities. But even 
in this case of high symmetry the questions left are still numerous as numerical 
simulations by M. Choptuik show [Q. He found very interesting properties, the 
so called echoing effect, for models which are in the parameter space of initial 
data near to the boundary which separates regular from singular spacetimes. 
In this paper conformal techniques are used to analyze the hyperboloidal initial 
value problem with scalar fields as matter models — for data near Minkowskian 
data the future of the initial value surface possesses a smooth future null infinity 
and a regular timelike infinity, for large data a smooth future null infinity exists for 
at least some time. In the second part of the introduction more about conformal 
techniques and their application for a mathematical description of asymptotically 
fiat spacetimes will be said. 

Although the primarily treated matter model is that of the conformally invariant 
scalar field, whose equations can be written as 

□0-:^0 = O (la) 

(1 - i^') = ((V,0)(V,0) - ^ V„V,0 - ^^,,(V'0)(V,0)) , (lb) 

the results obtained apply to a larger class of scalar field models, given by the class 
of actions (p2D, including the massless Klein-Gordon field, as shown in section ^ 
Note that an arbitrary factor can be absorbed into (p which changes the coefficients 
in ( [TBI) . My notational conventions are described in the appendix, the marks 
quantities in the physical spacetime (see definition |l|). The energy- momentum 
tensor for the conformally invariant scalar field can be written as 

fab = ( Va0) ( V,0) - ^ VaVt4> + \ 4>^ Kb ' \ Lb (( V'^) ( V,0) + ^ ^'i?) . ( Ic) 

The analytic investigation presented in this paper show the well posedness of the 
initial value problem in unphysical spacetime, which is a technical construct to 
"compactified" asymptotically fiat spacetimes in analogy to the compactification 
of the plane of complex numbers (R^) into the Riemann sphere (S^). 
One goal of this work was making myself familiar with the system in unphysi- 
cal spacetime as a preparation for numerical work showing that the conformal 
techniques are well suited for a numerical investigation of global spacetime struc- 



ture and gravitational radiation |21, 22|. To lower the computational resources 



required these calculations have been done for spherical symmetry. It is well 
known that spherically symmetric, uncharged vacuum models are Schwarzschild. 
The inclusion of matter removes that obstacle, the spacetime may evolve dynam- 
ically. Furthermore there is no gravitational radiation in spherically symmetric 
models. Therefore the matter model should also be a model for radiation. Scalar 
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fields with wave equations are choices for matter which model also radiation. The 
conformally invariant scalar field has been chosen since the matter equations are 
form invariant under rescalings of the metric and an appropriate transformation 
of the scalar field as the name already suggests. 

The scalar fields are interesting from the analytic viewpoint since for the first time 
conformal techniques could be used for matter models whose energy-momentum 
tensor has non-vanishing trace. 

1.2 Asymptotically flat spacetimes 

In this paper a geometrical, coordinate independent definition of asymptotical 
flatness along the lines suggested by R. Penrose will be used. A more thorough 
discussion of the ideas and the interpretation can be found at various places in 
the literature, e.g. [|r^, The definitions of asymptotical fiatness given in the 
literature differ slightly. The following will be used here: 

Definition 1 A spacetime (M, g^f^) is called asymptotically fiat if there is an- 
other "unphysical" spacetime (M, Qab) with boundary J and a smooth embedding 
by which M can be identified with M — J such that: 

1. There is a smooth function Q on M with 

^ Ia7> 9ab \m= ^^9ah- 

2. On J 

= and ^ 0. 

3. Each null geodesic in {M,gab) acquires a past and a future endpoint on J . 

Because of item ^ null geodesically incomplete spacetimes like Schwarzschild are 
not asymptotically flat. The next deflnition includes those spacetimes which have 
only an asymptotically flat part: 

Definition 2 A spacetime is called weakly asymptotically fiat if definition |i| 
with the exception of item |^ is fulfilled. 

Deflnition |l] and ^ classify spacetimes, they do not require that Einstein's equa- 
tion is fulflUed. One would like to know: 

Are they compatible with the Einstein equation with sources? Neither deflni- 
tion in nor ^ is in an initial value problem form: A given spacetime is or is not 
classifled as asymptotically flat. But for a physical problem one would like to 
give "asymptotically flat data" and have guaranteed that they evolve into an at 
least weakly asymptotically flat spacetime. 

Nevertheless the geometrically description was extremely helpful in analyzing 
asymptotically flat spacetimes and it can be successfully used as guideline to 
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construct a formalism which is better suited for analyzing initial value prob- 
lems. This method has been developed and applied to various matter sources by 
H. Friedrich ||lT], Il2|, O, |1^ Il5|, O] . In this paper it will be applied to general 



relativistic scalar field models. 

The idea is to choose a spacelike initial value surface in the unphysical spacetime 
{M,gab) and to evolve it. The problems to be faced are: 

For Minkowski space the unphysical spacetime (M, gab) can be smoothly extended 
with three points, future (i"*") and past timelike infinity, the end respectively 
the starting point of all timelike geodesies of (M, ^^j), and spacelike infinity 
the end point of all spacelike geodesies of (M, The point i° divides J' into 
two disjunct parts, future {J^~^) and past {J^~) null infinity. It is well known 
and has been discussed elsewhere that there are unsolved problems in smoothly 
extending a "normal" Cauchy hypersurface of M to if the spacetime has non- 
vanishing ADM mass. Certain curvature quantities blow up at reflecting the 
non-invariance of the mass under rescalings. 

By choosing a spacelike (with respect to Qab) hypersurface 5* not intersecting i° 
but J'^ (i7 ) we avoid the problems with i^. S is called a hyperboloidal hy- 
persurface — the corresponding initial value problem is called a hyperboloidal 
initial value (a detailed definition for the scalar field models is given in section ^ 
definition The domain of dependence D{S) of S will not contain the whole 
spacetime. The interior of S corresponds to an everywhere spacelike hypersurface 
in the physical spacetime which approaches a null hypersurface asymptotically. 
If is a light cone L then the domain of dependence of 5" is L. Therefore the 
hyperboloidal initial value problem is well suited to describe the future (past) of 
data on the spacelike hypersurface 5", e. g. a stellar object and the gravitational 
radiation caused by its time evolution. It is not well suited to investigate the 
structure near i^. 

But even for the hyperboloidal initial value problem there are "regularity" prob- 
lems at JT: Transforming the Einstein equation from physical to unphysical space- 
time an equation "singular" for Q = results. That problem is solved in this 
paper in analogy to H. Friedrich's work. A new set of equations for the unphysical 
spacetime will be derived, its equivalence to the Einstein equation on M proven. 
This new set of equations is used to prove the consistency of the hyperboloidal 
initial value problem for scalar fields with (weakly) asymptotical flatness and the 
existence of a regular future (past) timelike infinity for data sufficiently close to 
data for Minkowski spacetime. 



2 Regularizing the unphysical field equations 

A first attempt for equations determining (M, Qab) is the rescaled form of the 
field equation in physical spacetime. A closer look at the transformation of the 
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Einstein tensor under rescalings gab = ^'^9abi 

Gab = Gab + 2 {VaVb^ ' {V'V.Q) gab) + 3 ^]-2(V^^]) (V,^]) gab, (2) 

shows that this first attempt fails. Either there are terms proportional to 
and which need special care on the set X of points where f2 = 0, including 
J', which is part of M. Or alternatively, the highest (second order) derivatives 
of the metric, hidden in the Einstein tensor, are multiplied by a factor of fl'^ and 
then the principal part of the second order equation for the metric components 
vanishes on I. This behaviour of an equation will be called singular on I. 
In this section a system of equations without the singularity on I will be derived 
from the rescaled Einstein equation by introducing new variables and equations. 
The set of equations together with the equations for the matter variables may be a 
system with a very complicated principal part — as it is the case for a conformally 
invariant scalar field as matter model. A procedure is carried out to simplify the 
principal part to a form in which no equation contains both derivatives of matter 
variables as well as derivatives of geometry variables and the principal part of 
the subsystem for the geometry variables is the same as for the vacuum case 
("standard form"). All variables already present in the vacuum case are called 
geometry variables. 

It is shown that the procedure works for the conformally invariant scalar field. 
The procedure described does not use very restrictive assumptions — it is very 
general — and may work for most matter models, for which the unphysical matter 
equations can be regularized on X. 



2.1 The geometry part of the system 

According to the definition of asymptotical flatness (definition [^) the unphysical 
spacetime is connected with the physical spacetime through the rescaling 

gab\M=^^9ab- (3) 

This rule also determines the transformation of the connection and the curvature. 
Additionally the transformation of the matter variables $ under rescaling must 
be specified, 

$|^,= $[^,„f],l>]. (4) 

It is assumed that $ has a smooth limit on J , the rescaled equations for the mat- 
ter variables are regular on and there exists a tensor Tab-, which is independent 
of derivatives of Vt and derivatives of curvature terms, fulfills 

Tab \m= ^ ^ Tab, (5) 

^In the general case it is not known how to achieve that. 
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and has a limit on J 0. The conditions required may seem very restrictive but 
they can be fulfilled for Yang-Mills fields 0] and for the conformally invariant 
scalar field. 

The Riemann tensor will be split into its irreducible parts, the conformal Weyl 
tensor 

Cahc ='■ ^dab/, (6) 
the trace free part Rab of the Ricci tensor Rab and the Ricci scalar R: 

Rabcd = ^ dabcd + gc[aRb]d — 9d[aRb\c + -^gc[agh\dR, (7) 

A is used as an indication for trace free parts of 
tensors. 

The irreducible decomposition of the energy-momentum tensor is 

Tab = Tafe + ^ Qab T. (8) 

The irreducible parts transform under rescalings according to 

Tab = ^ '^Tab 

and 

T = n-^f. 

The vanishing of the divergence of Tab becomes 

= v^fab = fi^ v'^fab + VbT + n^T Vb^. (9) 

For energy-momentum tensors with non-vanishing trace equation as an equa- 
tion for the components of the irreducible parts of the energy-momentum tensor 
Tab is singular on X. Since should be in same way part of the matter equations 
problems in regularizing the matter equations are to be expected. 



2.1.1 A regular system 

The part of (H) proportional to is a pure trace, thus the singularity is 
absent in the trace free equation. A decomposition into the trace and the trace 
free part moves the worst term into one equation. 
From the rescaling rule for the Ricci scalar and tensor, 

R = n'^R + Qnv''Va^-i2{vn){Va^), (10) 

^From the definition of asymptotical flatness and the Einstein equation it follows that 
fl^^Tab has a limit on Q]. The faster fall off and the requirements on the form of Tab 
have technical reasons. 
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and 



Rab '■— Rab — -^QabR 

= Rab + 2n-^Vayb^-\^-\y'^Sl)9ab, (11) 

Gab = Tab, G = —R, and T = G it follows 

nR + Q vVa^ - 12 n-^ (v"fi) ( v„fi) = -n^T (12) 

and 

nRab + 2 VaVbft - ^ (Weft) gab = ft^ Tab- (13) 

Equation ( [T^ ) can be dealt with by the following lemma: 

Lemma 1 From R + f = (= ^) at one point, Gab = Tab (= and 
V Tab = R + T = follows everywhere. 

Proof: 

V'^fab = V'^fab + ^Vbf = 0. 

Combined with 

= V'Gab 

= ^"dab + \VbG 



Vfe(T + R)=0, 



gives 

i.e. T + i? is constant. 

Equation (|12]) will not be used any longer since V Tab = can be derived from 
the remaining equations, contract ( pOg] ) or see the discussion following (pl|). 
In the following the Ricci scalar R will be regarded as an arbitrary, given function. 
It fixes part of the gauge freedom on the transition from the physical to the 



unphysical spacetime as follows: The equations (|T2D and ([T3|) are invariant under 
rescalings {gab, ^) ^ {(jab, ^) := {Q'^gab, © ^) with 9 > 0. All the unphysical 
spacetimes (M, Q^gab, Q^) belong to the same physical spacetime (M, ^^^). 
Under the rescahng gab = gab, R and R are connected by 

6 VVa© = ei? - e^i?, (14) 

which is equation ([To|) where the covariant derivatives Va now corresponds to 
the unsealed metric. Solving ( p^ for a spacetime (M, gab) and data for 9 and 
G on a spacelike surface S we get at least locally a unphysical space time with 
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arbitrary Ricci scalar R. 

There is still conformal gauge freedom left as every rescaling with > and 

v"v„e = ^ei?(i-e2) (15) 

leaves the Ricci scalar unchanged. 

Equation (|T3D serves as regular equation for Vt. Substituting = j V^V yields 

VaVh^ = -^^Rab + i^gab + ^^^'fab, (16) 

which is a second order equation for f2. 

The next step is to find equations for the metric and the quantities derived there- 
from. Expressing the once contracted, second Bianchi identity (V[a-Rbc]d" = 0) in 
terms of Rab and dabc'^ results in 

y[aRb]c = - ^ {y[aR)9b]c - (Vdl)) dabc'' - ^ V ddab/ ■ (17) 

The once contracted second Bianchi identity in the physical spacetime, 

^dCabc"^ = -V[a(-Rfe]c - Q9b]cR)^ 

together with 



■'abc ) 



^VdCab/ — Vd(r2 ^Ca 

and the Einstein equation in physical spacetime provide us with an equation for 

rl 

^abc • 

"^ddab/ = V[aTfc]c — 3 (V[afi) T'tjc + gc[aTb]d ( Vfi) 

+ ^{Vla^)Tg,],+ ^n{VlaT)gb]c =: tabc- (18) 

We can now derive the missing equation for uj from the integrability condition 
for ( p!6D and by substituting (0): 

VaOO = -]^Rab V'n-^R Va^ - V.i? + fab 

-ln^Van)T-^n^VaT. (19) 

In the following in addition to the abstract indices (small Latin letters) frame 
(underlined indices) and coordinate indices (Greek letters) are used. The used 
conventions are explained in the appendix in more detail. 

Using Qa '■= Vaf^, the frame Cj", and the Ricci rotation coefficients 7"^- as further 
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variables, we get the following first order system of tensor equations for f2, 
u;, Ci", , Rab, and dab/- 

Va^ -^a = (20a) 

Vafif, + ^nkab - OOgab ' ^f^'fab = (20b) 
VatO + ^Rab^" + ^RQa + j^^^ aR ' l^^fab^" 

+^Q'QaT + ^Q'VaT = (20c) 
= (20d) 

-Rdiff ab/ — Ra\g ab/ = (20e) 
V[ai?fe]c + ^ {'^[aR)gb\c + ^^d^fe/ + tabc = (20f ) 

Vdrfafec'' - tabc = (20g) 

where ( pOd|) means vanishing torsion T"f,c, expressed in frame index form, 

= (ej(e/) - efe(e/)) e\ + Y'lk - 1% 



J\fr, 


— Efl a — 


■f^DSiab = 


E-Diiab = 


J^ujab 


= E^a = 


■^e°'bc = 


771 a 

- J^e be — 


\r d _ 

•' ^ 7 abc 


rp d 

-'-'7 abc 




E-Rabc = 




Ed abc = 



and ( P0e|) means that the curvature tensor in terms of the Ricci rotation coeffi- 
cients, in frame index form 

-Rdiff ijfe- = ej(7-ifc) - eii^j-jk) - 'ihjnl—jk + l-iml—ik 

x^.'rn. i I m I I rpm 

~r / vjl mk ~r / jjyf mk i mk-'- jii 

should equal the combination 

^ dabcd + gc[aRb]d — gd[aRb]c + -gc[agb]dR ='■ Ralgabc'^, 

which is the irreducible decomposition of a tensor with the symmetry of the 
Riemann tensor (|^) . Hence (|20d| ) and ( |20e| ) ensure that i?aig abc^ is the curvature 
tensor corresponding to the connection given by the Ricci rotation coefficients 
which again is the torsion free connection coming from the metric (frame). 



2.1.2 Complications by the matter terms 

The final goal is to use the terms Va^, Va^b, Va^^, {ej{ei/') — efc(ej'^)j e-^, 

ejil-ik) - eiiY-jj:), V[aRb]c, and V ddabc'^ in (|0) as principal part for the geom- 
etry variables of the system. I will call these terms left side of the equations, the 

^ The symbol E stands for equation, the first index reminds to the quantity for which a 
temporary equation wiU be formed by setting the tensor E equal to 0. The tensors providing 
the eventually used forms of the equations are named with M standing for null quantities. 
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remaining terms right side. The left side does not contain the complete principle 
part of the system yet as the energy momentum tensor Tab and its derivatives 
V[aTf,]c may contain derivatives of the matter and geometry variables. 
In the case of the conformally invariant scalar field the field equation (|l|) remains 
invariant under the rescaling 

(f) = 

i.e. 

□0_ = 0. 

The physical energy-momentum tensor Tab fulfills the assumed properties, 



Tab — 



(Va0)(V60) - ^0VaVb0 + ^(f)^Rab - \9ab ((V^0)(Vc0) + R 



n^T 



ab- 



The mentioned complications in pO|) by the right sides are now obvious. 
Firstly V[aTfc]c contains V[a'Vb]Vc4> terms which are eliminated with the iden- 
tity V[aVb]Vc0 = |-Rafec'^Vrf(/). To get rid of the second and first order derivatives 
of (j) we use the first order system 



^f^a = '^a(f>-4>a = (21a) 

1 

— ( 

4 



A/'o./.afe = Va0fe - 0afe - 70/5'a6 = (21b) 



^f^^=K-^^ = (21c) 


1 1 

A/'dD^ abc = V[a06]c + q ^ [^R + R 4>[a)gb]c - -Ralg abc'^^d = (21d) 
^fD^^a = Va<Pb'-]:{^'^aR + R4>a)=0. (21e) 

o 

for the variables (pa, the trace free symmetric tensor (pab and the trace (pa"'- The 
system is derived from Va(n0— -f^) = 0. System (|2T|) also serves as matter 
part of the system for the unphysical spacetime. tabc is now written in a form 
which does not contain any derivatives of matter variables explicitly. 
V[aTb]c and thus tabc still contain derivatives V[aRb]c of the trace free Ricci tensor. 
By combining (|2m]) and (|20g|) the derivatives of Rab and dabc'^ can be decoupled. 



0|) and ( |20g| ) become 



E'nabc = V[a-R6]c + ^ ( V[ai?)fi'6]c " ^ddabc'^ + ^^labc = (22) 

and 

E'dabc = ^d.dabc'^ — "^abc = 0, (23) 
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with 



triabc 



1 - i fi202 



9c[a 



— (f)[agb]cR + ^ ^ (f)^dabc'^^c 



1 

4 

l2j d( 



+ -(t) Qdab/(j)d + ^4> gc[aRbf<Pd - 



-3% 



^b]<Pd 



+ ^ (p^Rb]c + ^ 4>'^gb]cR - ^ s-feic 4>'^4>d 
id + ^ <P'^Rb]d 



'[aRb]c 
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Note that mahc may become singular for l — ^fl'^(f)'^ = l — ^(j)'' = 0. In the Einstein 

equations for the physical spacetime (|r^) Rab carries a factor 1 — ^0 too. We 
will need later that 

1 2 f 2 d \ 

m abc '■= tabc — ^abc = —^^<P ( A/'r afcc + g ^ "^[ald] fi'6]c j , 

where A/'aafec is the null quantity representing the final form of the equation 
for Rab (p5|). The final form of the equation for dabc"^ is obtained from (^) by 
replacing E'd abc = with 



A/'d afec := E'dabc + ■^fn[a\d\''' gb]c = 0. 



(24) 



This gives Afd abc the same index symmetry properties as the Weyl tensor. That 
replacement does not change the equation since mab'' = as will be seen later. 
Analogously we replace (^2]) with 



■f^Kabc '■— E'-Elabc — - ^^[a\dfgb]c 



0, 



(25) 



the contraction M-Rab^ = is then the contracted second Bianchi identity. 



3 Evolution equations and constraints 



In the following I will assume a system A/" = of the form (|20a|) - ( pOe| ), (p^ 



(P^D, the geometry part, and a matter part, in the case of the scalar field model 



system (|21|) . mabc and tabc are assumed to differ only by terms expressible as null 
quantities. The energy-momentum tensor Tab, its derivatives V aTbc, ^abc, and 
tabc are assumed to be expressed in variables and thus do not contain any explicit 
derivative of variables. By these assumptions the principal part of the system 
has block form, the geometry block and the matter block. The two blocks are 
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coupled through the right sides. 

In this chapter the system A/" = will be reduced to a system of symmetric 
hyperbolic time evolution equations, the subsidiary system. Sufficient conditions 
for the equivalence of the subsidiary system and Af = are given as conditions 
on rriabc- If the system can be put into the described block form there do not 
arise any more conditions from the geometry part of the system for any matter. 
The explicit carry out is technical and lengthy, the idea can be summarized as 
follows: All the equations of the system are regarded as null quantities. By requir- 
ing the vanishing of some of these null quantities and by choosing an appropriate 
gauge condition for the coordinates and the frame we get a symmetric hyperbolic 
subsidiary system of evolution equations. Which null quantities to choose can 
best be seen by a decomposition into the irreducible parts in the spinor calculus 
as performed in [ 1^ . 



The solution of this symmetric hyperbolic subsystem exists and is unique. To 
complete the proof we must show that the solution obtained in this way is consis- 
tent with the rest of the equations, i.e. that all null quantities remain zero if they 
are initially zero ("propagation of the constraints"). For that purpose a symmet- 
ric hyperbolic system of time evolution equations for the remaining null quantities 
is derived. Sufficient conditions for the propagation of the constraints are ffistly 
the homogeneity of the evolution equations for the remaining null quantities in 
the null quantities since then the unique solution of these evolution equations is 
the vanishing of all null quantities for all times if they vanish on the initial surface 
and secondly that the domain of dependence of 5" with respect to the equations 
for the propagation of the constraints is a superset of the domain of dependence 
of S with respect to the subsidiary system. 

3.1 A symmetric hyperbolic subsystem of evolution equa- 
tions 

Introducing a timelike vector field not necessarily hypersurface orthogonal and 
its orthogonal projection tensor hat '■= Qab — tatb/(tct^) allows to split the system 
of equations into two categories, the equations containing time derivatives and 
those containing no time derivatives (the constraints). The equations with time 
derivatives provide a under /overdetermined system of evolution equations. 
The system is overdetermined since for some quantities there are too many time 
evolution equations, e.g. there are 12 time evolution equations from M-Rabc = 
for 9 independent tensor components. 3 equations are a linear combination of 
the other 9 equations and the constraints. An irreducible decomposition of the 
tensors A/" is a systematic way to analyze these dependencies. Since all the 
types of tensor index symmetries appearing in the system have been thoroughly 



investigated in I will only state which combinations are needed. 



The system is underdetermined since there are 10 time evolution equations for 
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the frame and the Ricci rotation coefficients missing. By adding 



and 



eh g'-^ e/ ( VaCfc') = -F, = 7^^^ (26) 



dkfM + yMf, + 7^ 7/fc^ - 7^- Ilk- = F% (27) 

the system becomes complete. The freedom of giving ten functions corresponds 
to the freedom of giving the lapse and the shift to determine the coordinates and 
the six parameters of the Lorentz group to determine the frame on every point. 



The gauge freedom is discussed in full detail in |T^, [T^ • 

A choice which makes the system especially simple for analytic considerations is 
a Gaussian coordinate and frame system defined as follows. Give on the spacelike 
initial value surface S coordinates x^,^ = 1..3, and 3 orthonormal vector fields 
Gi'^i = 1--3 in this hypersurface. The affine parameter of the geodesies of the 
hypersurface orthonormal, timelike vector field eo° defines the time coordinate 
= t. The spacelike coordinates are transported by these geodesies into a 
neighbourhood of the initial surface. By geodesic transport of Cj", i = 1 . . . 3, and 
Co" a frame is obtained in this neighbourhood. By construction we have 

eo° = 1, e/ = for = 1..3 

and 

Y-ok = 0. 

It is well known that Gaussian coordinates develop caustics if the energy momen- 
tum tensor fulfills certain energy conditions, see e.g. |2^, lemma 9.2.1]. In the un- 
physical spacetime the fl terms provide a kind of unphysical energy-momentum 
tensor. Whether this energy-momentum tensor fulfills the energy-momentum 
conditions is a difficult question and not known to the author. Nevertheless the 



coordinates develop caustics as has been shown by numerical calculations P2 |. 
The following combinations give a symmetric hyperbolic system for the remaining 
variables as can be deduced from the considerations in 11121 : 



Arno = 0, (28a 

A^Dnof, = 0, (28b 

AT^o = 0, (28c 

^f.^o = 0, (28d 

AT^oi/ = 0, (28e 

^^^'A^Kiafe = 0, 1=1,2,3, (28f 

A/'rm = 0. 1 = 1,2,3, (28g 

^fKOiJ_ + ^f■ROli = 0, (i,j) = (1,2), (1,3), (2,3), (28h 

AAd 212 - ATd 313 +Ard 202 -ATd 303 = 0, (28i 
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— A/ H 1 09 + A/ H 1 91 


= 0, 




(28i) 


A/'d 101 = 0, 






(28k) 


-A/ H 1 09 -A/ H 1 91 — 


0, 




(281) 


— TV H 91 9 ~l" -A/ H '^1 3 


~l" -A/ H 909 — -/V H 303 


= 0, 


(28m) 


* (J. Z_LO 1 * CJ. O-LZ 1 


yV H 90*^ H~ yV H '^09 = 


0, 


(28n) 


— A/ H 1 OS + A/ d 1 


= 0, 




(28o) 

V / 


-A/'d 123 = 0, 






(28p) 


— A/'d 103 — A/'d 131 


= 0, 




(28q) 


A/ H 91 ~l~ A/ H 31 9 — 




0, 


(28r) 


AT^o = 0, 






(28s) 


A/" Dfliofe = 0, 






(28t) 




1= 1,2,3, 




(28u) 


A/'dD()!>0m = 0, 


1 = 1,2,3, 




(28v) 


A/" DDflioij' + A/'ddc^ 


oji = o, = 


= (1,2), (1,3), (2, 3), 


(28w) 


A/'Dn(/)0 = 0. 






(28x) 



To see that the system is really symmetric hyperbolic one has to write down the 
system explicitly. By an appropriate, in the explicit form of the system obvious 
definition of new variables, the system has the structure 

A^dJ + j2A^d,.l + kil, x'^) = 0, (29) 

i=l 

with a diagonal matrix At, which is positive definite for 1 — lfi^0^ > 0, and sym- 
metric matrices A^i. f is the vector build from the variables. All the remaining 
equations are linear combinations of (pSj) and constraints. Since the explicit form 
of the constraints is not needed I do not list them. 

As the entries in At coming from Afn = and TVd = vanish for 1 — = 

the following results apply only if < 4 everywhere on the initial value sur- 
face S and thus in a neighbourhood of S. The physical Einstein equations have 
a corresponding singularity for 1 — ^0 =0 (see equation p!E|) . 

3.2 A sufficient condition for the propagation of the con- 
straints 

According to the analysis in Jl^, involving the left hand side of the following 
identities, a symmetric hyperbolic system of evolution equations for the remaining 
null quantities can be extracted from: 

"^[aAfn b] = -^T^afe V,fi - ATjin [ab] (30a) 
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[a-'V on b]c — 

1 1^ 1 1- 

[<^9b\c + ^f^^"^[a|d|'^fl'fe]c + ^^Wm [a|d|'^fi'fe]c (30b) 



V[aA/'a,6] = 

+hl Q[6m„]/ + ^fi nW^abc (30c) 



2"^ 7 [abc/Rfd + 2"^ 7 [a&|d]'^-^c]/ " [ab^ \f\Rc]d " ^7'-''[a&( V|/|-R)5'c](i 

+A/^Dn [a\f\dbc]d'^ + ^-^{Md [ca\d\9b]f — J^d [ca\f\9b]d) + J^Si [a"^6c]d 

2 1 - 1 - 1 

--^J^Silarriblf/gcjd - -^^"^J^-ylabc/Tfa - -Vt'^M^[ab\d\^Tc]f - -Vt^T^ab^ \f\T c]d 

.1 1 

-SQjV^DSl [abTc]d + f^A/'on [a 7'c|/|fl'6]d + -^^A/'dO [ab9c\d + J^^A/^n [a(VbT)g'c](i 

m [ca 9b]d 

2 

-QV[aA^m6c]d + 2n[amc\f/ 9b]d " g ^^( V[amb|/|-'')5fc]d (30d) 



VWda6c = 

1111 

ab} + 2"^ 7 d/C^abc 

-^T^^'^V/46cd - V^maftc + \y[bma]c'' (30e) 

V[aA/'e '^ftc] = J^-r [abcf + -^e \ab^ e ''c]/ (30f ) 

^ [fJ^ -y ab]cd = 

-T'^l/a-Rdiff 6]gcd - A/'n [fdab]cd " ^^(A/'d [6/|c|5'a]d " d [bf\d\9a]c) 

—J^n [fb\d\9a]c + J^R [fb\c\9a]d (30g) 
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V[aAr^ b] = -^r",feVc0 - ^fu^ [ab] (30h) 



1 d 1 d 1 

VfaA^D^feJc = --^-/ahc (l>d — -T ab^ d't'c — -^HTlipabc — --^ Dn4,[a9b]c (30i) 

Va^f^^ = ATon^a - ^I^^a (30j) 



V[aA/'DD(/)6c]d = 

1 f-^ 1 r If " 1 

-jN'^abJ (t)fd + 2-^7[af'M| V'c]/ - 2^ + g A/",^ [a( V6-R)5'c]d 

If 1 1 r 1 . 

— ^<pT'[ahgc]d^fR + -^Rf^T)<P[ab9c]d+ -^T^ [abRdiS c]gd <P f — ^R[bc\d\ J^r>(f)a]f 

+ li'^la^f^bc]/)(pf, (30k) 

and 

V [a A^D □</.{)] = —-T'^ab^ d4>c — -A^(^[aVb]-R — —(jyT^ab^cR " --RA^D,/, [afe] • (301) 

The last term in ( p01<| ) is homogeneous in null quantities as can be seen from 
The deviation of these equalities is even more lengthy than the equalities 



itself, but the essential ideas behind it can already be seen in the deviation of the 
first: 

V[aA/'n6] = ^[aSb]^ — V[„fi6] 

= —-T^'ab^ Sf' — A/'dH [ab\ + -R[ab\^ " ~ ^dlab], 

with the last three terms vanishing since the tensors are symmetric. Note that 
vanishing of the torsion, T^ab = 0, and 2 V[aVb]^^c = Rabc^^d cannot be used since 
they only hold if both the time evolution and the constraint equations for the 
frame Cj^ and the Ricci rotation coefficients 7"^- hold everywhere. 
A sufficient set of conditions for homogeneity of the system derived in the null 
quantities is 

^ab^ = mod A^ , (31a) 
2 

V'mabc + - V[am6]/ = mod N , (31b) 
V[amfe]/ = mod U , (31c) 

and 

^^V[aA/"m6c]<i = /V[aArR6,]rfmodAr, / ^ -1. (31d) 
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A straightforward but long calculation shows that these conditions are fulfilled by 
the conformally invariant scalar field with / = — Equation (|30d| ) becomes 

singular for 1 — = 0. 

The very technical integrability conditions (|3T|) have a very simply interpretation. 

Replacing niabc with tabc — they only differ by null quantities — the conditions 

~ b ~ ~ ~ ~ 

(|31aH31cj ) reduce to V Tab = and V^V Tac = 0. Condition ( |31d|) is only of 

technical nature, it gives the principal part a simple block form. 

From the considerations in |0 also follows that the domain of dependence of 5* 

with respect to the evolution equation of the constraints includes the domain of 

dependence of S with respect to the subsidiary system. 



4 The hyperboloidal initial value problem 

So far a system of equations {M = 0) has been derived which contains for at least 
one choice of gauge a symmetric hyperbolic subsystem of evolution equations. 
The remaining equations in the system — either constraints or a combination 
of constraints and time evolution equations — will be satisfied for a solution of 
the evolution equations, if the constraints are satisfied by the initial data. If 
both, the time evolution and the constraints, are fulfilled, (M, g^h, (p) is a weakly 
asymptotically fiat solution of the Einstein equation. This follows from the way 
the system A/" = for the unphysical spacetime has been derived. 
The essential points in the proofs of the theorems in chapter 10] are the sym- 
metric hyperbolicity of the subsidiary system and the form (|20bD of the equations 
for Q. Therefore the same techniques can be used and the proofs will not be re- 
peated. The difference to the model treated here lies in the derivation of the 
subsidiary system and the proof of the propagation of the constraints, which has 
been done in the previous chapters. 

4.1 The initial value problem 

We consider the following initial value problem: 

Definition 3 A "hyperboloidal initial data set for the conformally in- 
variant scalar field" consists of a pair {S,fo) such that: 

1. S = S {~\ dS is a smooth manifold with boundary OS difjeomorphic to the 
closed unit ball in M'^. As coordinates on S the pull backs of the natural 
coordinates on M'^ are used. 



2. fo is the vector f of functions in system written in the form ^2^) at 
initial time to- 
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3. The fields provided by fo have uniformly continuous derivatives with respect 
to the coordinates of S to all order^. 



4- On S: > 0. On dS: ^1 = and Vafi is a future directed null vector. 
5. The fields provided by fo satisfy the constraints following from Af = 



and \2j\) and the gauge conditions. 



A point which deserves special notice is the existence of a hyperboloidal initial 
data set. The proof that those data exist has to overcome two problems. 
Firstly, the regularity of the solution on dS which is the consistency of the data 
with asymptotical flatness. For scalar field data with compact support regularity 
conditions are given in [0, 0], which are sufficient for the existence of a solution 
of the constraints near dS. 

Secondly there is a problem with a possible singularity of equations in A/" = at 

1-1^2^2^0. 

L. Anderson and P. Chrusciel are preparing a paper analyzing both problems M. 



4.2 Theorems 

The "theorems" will be given in a form not containing every technical detail, since 
these technical details would make them lengthy and can be easily deduced from 
the theorems in |jT^ by replacing the Yang-Mills matter with the (conformally 
invariant) scalar field. 

Since the constraints of A/" = propagate we have: 

Theorem 2 Any (sufficiently smooth) solution of the subsidiary system satisfy- 
ing the constraints on a spacelike hypersurface S and 1 — | > defines in 
the domain of dependence with respect to gab of S a solution to the unphysical 
system. Thus (M, 0) is a weakly asymptotically flat solution of the Einstein 
equation. 

Since the evolution equations are symmetric hyperbolic a unique solution of the 
initial value problem exists for a finite time. From the combination with theorem 
(0) follows: 

Corollary 3 For every regular solution of the constraints on S with 
1 — \ VL^(t)'^ \s> exists locally a unique, weakly asymptotically flat solution of 
the Einstein equation. 

For the Minkowski space we can extent S and the solution of the constraints 
beyond dS to 5" and get a solution in the unphysical spacetime which extents 
beyond i^ . The continuous dependence of the solution of symmetric hyperbolic 

"^The assumption about the smoothness of the data can certainly be weakened from C°° to 
C" for sufficiently large n but then more technical effort would be needed in the proofs. 
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systems on the data and the form of ( |20a| ), ( pOb| ) and ( pOc| ) (see the proof of 



theorem (10.2) in [|T^) guarantees that there is a solution covering the whole 
domain of dependence of 5*. Furthermore the proof there shows that {p \^{p) = 
0} has an isolated critical point where all future directed timelike geodesies 
of (M, g^h) end, thus: 

Theorem 4 For a sufficient small deviation of the data from Minkowskian data 
the solution of theorem possesses a regular future null infinity and a regular 
future timelike infinity. 

5 The conformal equivalence of the scalar fields 

This section shortly reviews the equivalence transformation between spacetime 
models with scalar matter under the viewpoint of solving hyperboloidal initial 
value problems. Other aspects of this equivalence transformation, especially the 
generation of exact solutions, have been studied in [0, |], |, |2^, 

5.1 Local equivalence of solutions 

Spacetime models (M, gab, 4>) with scalar matter described by the action 

5 = / \A{4>)R - i?(0)(V„0)(V>)l (-^)^ d^x (32) 

will be considered. Boundary terms in the action have been omitted, g is the 
determinant of ^^y. 

By varying the action S with respect to and g^b the following field equations 
result: 

2 



{VaVbA{4>)) + [VV,A{4>))gab = 0. (33b) 



A and B are assumed to be C°° functions. For B ^ the principal part of ( p3a| ) 
does not vanish and thus ( |33a| ) is a wave equation. For that reason I assume 
5(0) > e > for every 0. (|33b| ) is a second order equation for the metric if 

In the spacetime region i7 := |s G M| sign(y4($)) > V$ G [0o,0(x)]| the 
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transformationQ 



3 fdA" 



^^2 

dab = A gab 



+ AB 



(34a) 
(34b) 



gives a solution of the system (p3| ) with a massless Klein- Gordon field as matter 
model corresponding to the choice {A, B) = (1, 1) and the equations 



□ = 
Rab - -^R9ab= Tab[4>] 



with energy momentum tensor 



fab[^] = (V.0) (V,0) - ^ (V,0) (V'0) g,,. 



(35a) 
(35b) 

(35c) 



From the assumptions about A and B follows that the corresponding Klein- 
Gordon field will be unbounded approaching the part of the boundary of H 
where ^4(0) — > 0. The singularity in the Klein-Gordon field shows up at least in 
a singularity of the equations for (M, gab, 0). 

For two of the scalar fields in the above class the field equations are very special, 



the already mention massless Klein-Gordon field ([35| ) and the conformally in 



variant scalar field 0, {A,B) = (1 — ^0 , |) (0 can be rescaled by an arbitrary 
factor). 

The first, because the set of equations in the physical spacetime becomes re- 
markable simple and has been analyzed intensely with analytical (e.g. [^) and 
numerical (e.g. 0|) methods for spacetimes with spherical symmetry. 
The second, yielding the equations (U), because the matter equations are invari- 
ant under rescalings gab = (jab = f2~^0. 
The transformation between the two special cases is 

= \/6arctanh^ (36a) 
1 -2 ^ 

~9ab = (1 - ^0 ) 9ab (36b) 

which is a bijective mapping from g] — 2, 2[ to g] — oo, oo[. 
Due to the following diagram, illustrating the above described relations, it is evi- 
dent that there is a variable transformation regularizing the unphysical equations 
for the Klein-Gordon field: 



^The choice of the parameter reflects gauge freedom. Models where there is no (f>Q with 
A{(j>Q) > will not be considered. 
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physical 


conformal field (f) 


CITlP PPflTTIP 
O LJ CX Cu 1 1 1 iC 




-2,2[ 






= ^/Vt 






gab = ^'^ 


unphysical 


conformal field 0, 


spacetime 


reg. equations 



= f{4>) 



9ab 



^^9ab 



9ab 



KG field 




oo, oo[ 


= l/n 








KG field 0, 


sing, equations 



By mapping an arbitrary scalar field with action (32) to the Klein-Gordon field 

and then to the conformally invariant scalar field regular equations for are 
obtained. 



5.2 The hyperboloidal initial value problem 

Since A{(j)) > on all scalar field models connected by transformation (|3^ to 
a hyperboloidal initial value problem with a conformally invariant scalar field as 
matter source are weakly asymptotically fiat. 

For a massless KG model there is a one parameter gauge freedom in the scalar 
field. If is a solution, then so is 0+0q with 0q = const, as the energy momentum 
tensor depends on derivatives of only. This can also be seen by mapping a Klein- 
Gordon model to a Klein-Gordon model with 0q 7^ and (|3^ . The analogue 
holds for every considered scalar field model. For the hyperboloidal initial value 
problem = Q(f), therefore vanishes at JT", fixing the gauge in 0. 
In definition (|) 1 - ifi202 |^> was assumed. But with the Bekenstein black 
hole a weakly asymptotically fiat solution is known where ^4(0) vanishes on 
a regular part of the spacetime. In this case the transformation gives a possible 
extension of a massless Klein-Gordon scalar field solution beyond a singularity - 
the Klein-Gordon field and the metric g^^j, degenerate there. It is a pleasure for 

me to thank Helmut Friedrich, Bernd Schmidt, and Jiirgen Ehlers for the very 
helpful discussions during the grow of this work which is part of my Ph. D. thesis. 
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A Notation 



The signature of the Lorentzian metric gab is (— , +, +, +). 

Whenever possible I use abstract indices as described in chapter 2]. Small 
Latin letters denote abstract indices, underlined small Latin letters are frame 
indices. For the components of a tensor with respect to coordinates small Greek 
letters are used. The frame (gfjr) is constructed from the coordinates x^, Cj" 
denotes an arbitrary frame. In this notation Va is a covector, Vi a scalar, namely 
Va e^. 

v{f) is defined to be the action of the vector on the function /, i.e. for every 
covariant derivative Va- t{f) = f^Vaf- 

The transformation between abstract, coordinate, and frame indices is done by 
contracting with e," and Cj'^. All indices may be raised and lowered with the 
metric qab and the inverse g"^^ . g^^ gcB = ^^b, A and B are arbitrary indices, 
e.g. Cia = gab^i and e\ = g'-^-eta- 

For a frame e^" and a covariant derivative Va the Ricci rotation coefficients are 
defined as 

:= e,^V,e/. 

From this definition follows 

e^eib (Vat") = e,(^^)+7^^^-. 
With respect to a coordinate frame e^" = (gfjr) the components •j^fj,^ are the 
Christoffel symbols F^^,^. 
The torsion T'^^c is defined by 

V.VJ-VfeVa/ = -T^.feVJ, 
the Riemann tensor Rabc'^ by 

Va Vfeti^c — VfeVal^c = Rabc^ — T'^ ab Vrfl^c- 

Contraction gives the Ricci tensor, 

-'^ab ^acb ; 

and the Ricci scalar 

R = Rab g"" ■ 

The Einstein tensor is given by 

Gab = Rab — 2 RQab- 



The speed of light c is set to 1 as the gravitational constant k in Gab = i^Tab- 
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